Theory of a quantum critical phenomenon in a topological insulator: (3+1) 
dimensional quantum electrodynamics in solids 
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We study theoretically the quantum critical phenomenon of the phase transition between the 
trivial insulator and the topological insulator in (3+l)-dimensions, which is described by a Dirac 
fermion coupled to the electromagnetic field. The renormalization group (RG) equations for the 
running coupling constant a, the speed of light c, and electron v are derived. The almost exact 
analytic solutions to these RG equations are obtained to reveal that (i) c and v approach to the 
common value with combination c 2 v being almost unrenormalized, (ii) the RG flow of a is the 
same as that of usual QED with c 3 being replaced by c 2 v, and (iii) there are two crossover mo- 
mentum/energy scales separating three regions of different scaling behaviors. The dielectric and 
magnetic susceptibilities, angle-resolved photoemission spectroscopy (ARPES), and the behavior of 
the gap are discussed from this viewpoint. 

PACS numbers: 73.43.Nq, 71.10.-w, 64.70.Tg 



I. INTRODUCTION 

In solids, the electronic states are described by the 
Bloch wavefunctions with the energy dispersion e n (k) 
where n being the band index and k the crystal mo- 
mentum. The velocity of electrons given by v n {k) — 
de n {k)/dk is usually much smaller than that of light c. 
Therefore, the Lorentz invariance is terribly broken and 
hence many of the beautiful results in quantum electro- 
dynamics (QED)i are not applicable to the Bloch elec- 
trons in solids. The smallness of the factor v n (k)/c <C 1 
naturally leads to the gauge choice (i.e., Coulomb gauge) 
where the scalar potential gives the Coulomb interaction 
without retardation while the electron-electron interac- 
tion through the transverse part of the vector potential 
A is often neglected. The latter is often treated as the 
external electromagnetic field for the probe of the elec- 
tromagnetic response of the system. This gauge choice 
is regarded as the "physical gauge." For example, one 
can discuss the physical meaning of the Green's func- 
tion G(k,uj) in this gauge where the quasi-particle cor- 
responds its pole structure. Angle-resolved photoemis- 
sion spectroscopy (ARPES) is also formulated in this 
gauge [i.e., ARPES intensity is proportional to the elec- 
tron spectrum function — ImG(fc, w)} 2 -. 

While the nonrelativistic quantum mechanics is basi- 
cally justified for the electrons in solids, there are some 
cases where the Dirac fermions appear in the electronic 
band structure. A representative case is graphene, a two- 
dimensional sheet of carbon network with hexagonal lat- 
tice, where the 2x2 Dirac spectrum near K and K' points 
describes the low-energy physics^. Another example is 
Bi, which is described by 4 x 4 Dirac fermions and shows 
the enhanced orbital diamagnetism^. Recent advances in 
this field are the discovery of the topological insulator and 
its associated quantum phase transition^. The relativis- 



tic spin-orbit interaction (SOI) rearrange the spin states 
to yield the "twist" of the Bloch wavefunctions in the 
first Brillouin zone. This twist is characterized by the Z% 
topological integer. In general, topological integers can 
change only discontinuously when the gap closes, which 
can be described by the local Hamiltonian in k space. 
When the inversion symmetry exists, the effective Hamil- 
tonian near this quantum phase transition is the Dirac 
Hamiltonian expanded around the time-reversal invari- 
ant momentum (TRIM) fco (fco is equivalent to — fco). In 
this case, the orbitals and spins are coupled to form the 
4x4 Dirac Hamiltonian and the sign change of the mass 
m corresponds to the quantum phase transition between 
trivial insulator and topological insulator. This story is 
actually realized in the materials such as BiT^Si-^Se^ 
by changing the concentration 

The effects of the electron-electron interaction on the 
Dirac electrons are also extensively studied 9-12 . For 
graphene, it has been revealed that the electron speed 
v is renormalized to increase logarithmically by the 
long-range Coulomb interaction, while the coupling a is 
marginally irrelevant 9 -. When the exchange of the trans- 
verse part of the vector potential is taken into account, 
the velocity v saturates to that of light c (i.e., the Lorentz 
invariance is recovered) and a remains finite in the in- 
frared limit. This leads to an intriguing non-Fermi liquid 
state in (2+1) dimensions 10 . For the (3+l)-dimensional 
[(3+l)D] case, the Coulomb interaction also gives the 
logarithmic enhancement of the velocity v and the cou- 
pling constant a is marginally irrelevan t 11 : 12 . The disor- 
der potential is irrelevant perturbatively, while the strong 
enough disorder drives the system toward the compress- 
ible diffusive metal (CDM)iiii2. However, the effect of 
the transverse part of the vector potential in (3+1) di- 
mensions has not yet been studied to the best of our 
knowledge. 

In this paper, we study the quantum critical phe- 



nomenon of topological phase transition in (3+1) dimen- 
sions. The Coulomb interaction as well as the transverse 
current-current interaction are considered. 



II. DIRAC FERMIONS IN (3+1) DIMENSIONS 
IN ABSENCE OF LORENTZ INVARIANCE 
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FIG. 1. Feynman diagrams considered here: (a) self-energy, 
(b) polarization, (c) vertex. 



We start with the following Lagrangiani^: 

C =-4>h°po -vf-p- m)tp + -(eE 2 - —B 2 ) 

2 [i 

— etpj Q ipA — e—'ipj a ipA a , 
c 



(1) 



where a is a spatial index (a — 1, 2, 3) and 7 Q p Q = —j-p- 
For the moment, we consider the critical point (i.e., 
to = 0). The renormalization of the mass m will be dis- 
cussed later. The speed of light in material c and in vac- 
uum Cvacuum = 3 x 10 8 m/s are related through the per- 
mittivity e and the permeability /i by c 2 = c 2 acuum /(e/Li). 
We use a (H ) metric. The electric field and mag- 
netic field are represented in terms of the photon field 
as 

E = ——-VA , B = -VxA. 
c at c 



The electron propagator Go(p), the photon propagator 



Dq"(p) and the vertex are given by 



G (p) = 



7°Po + vj a p a + i0 ' 



vertex = —ie~/° or — ie— 7°. 

c 



(2) 
(3) 
(4) 



Here we employ the Feynman gauge because physical 
quantities are independent of gauge choice. 



B. Perturbative renormalization group analysis 

Calculations are performed by using dimensional regu- 
larization not to violate the gauge invariance of the the- 
ory. We set the space-time dimension d = 4 — e to reg- 
ularize divergences. The self-energy £(p), polarization 
U^iq), and the vertex correction 5T fJ, (p',p) all to one- 
loop order (Fig. [T]) are obtained as follows: 
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Comparing the result of the vertex correction Eq. ([7]) 
with the self-energy Eq. ([5]), we can confirm that the 
Ward-Takahashi identity is satisfied to one-loop order. 

The diverging quantities appearing through the calcu- 
lation of the one-loop diagrams are absorbed by rescal- 
ing some quantities. We can write the renormalized La- 
grangian in the form 

C =iiZ 2a °k + Z 2s v<y a p a )i> + ^{Z 3e eE 2 - Z 3m ±-B 2 ) 



- eZ lt ip^°ipA - eZ ls -ipj a ipA a . 
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Then we obtain the following RG equations using a mo- 
mentum scale k: 
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(11) 



The coupling constant a is defined by a = e 2 / (ec 2 v). The 
numerical solutions for the RG equations are shown in 
Fig. [2] for the initial (bare) values of vo = 0.01, Co = 0.5, 
and a a = 1. 

The result shows some important features. First, we 
can see that the quantity c 2 v is almost constant for all 
momentum scales and remains CqVq. This fact enables the 
approximate but accurate analysis of the scaling func- 
tions as described below. Second, the speed of elec- 
tron v and that of photon c approach to the same value 
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FIG. 2. (Color online) Solution to the RG equations for v, 
c, and q. We set the initial values vo — 0.01, co = 0.5, 
and ao = 1. For v, c, and a, the analytic solutions (solid 
lines) show a very good agreement with the numerical solu- 
tions (points). The dashed line for c 2 v, obtained numerically, 
is almost constant for all momentum scale. 



Coo = (co^o) 1 ^ 3 i n the infrared (IR) limit. Third, the cou- 
pling constant a becomes small in the IR region, which 
justifies our perturbative RG analysis. Therefore, the 
quantum critical phenomenon of 3D topological insula- 
tor is an ideal laboratory to study the QED in a solid, 
even though the Lorentz invariance is broken to a large 
extent in the original (bare) Lagrangian. 



C. Analytic solutions 

Now we study the solution to the RG equation in more 
detail. The approximate relation c 2 v = CqVq makes the 
analysis much easier, and we can obtain the analytic so- 
lution. By replacing c 2 v by CqVq, the RG equation for 
e 2 /e [Eq. (jlip] is exactly the same as in the conventional 
QED. Therefore, the RG equation for the coupling con- 
stant a is 



da 

d,K 



a 

6^' 



and its solution is obtained as 
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(13) 



where kq is the momentum cutoff. This approximate 
solution fits the numerical solution very well as shown in 
Fig. [2] The precision of the analytic solution is discussed 
in Appendix. 

With (x(k) being obtained, one can solve the RG equa- 
tion (fTOl) for c as 
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FIG. 3. (Color online) Analytic solutions to the RG equations 
for the permittivity e and the permeability fi. The character- 
istic momentum scales are different for e and fi. 



and v(k) = CqV /[c(k)] 2 . These analytic solutions are 
again compared with the numerical solutions in Fig. [21 
and a good agreement is obtained. 

Here we can identify the two momentum scales, Ki and 
K2- Ki is defined as the scale where the renormalization 
effect becomes appreciable, [i.e., ^ln(^) = 1]. The 
second one k,2 is defined as c{K2) — v(k2) (i.e., the two 
velocities approaches to each other). These two scales 
are estimated as 
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and K2 -C Ki <C kq, assuming ao/(67r 2 ) 1 and vq/cq <C 
1. These two momenta separate the three regions: (i) 
perturbative region Ki <C k <C kq, the renormalization 
effect is small and perturbative; (ii) non-relativistic scal- 
ing region k<i <C k <C «i, the renormalization effect is 
large, while c(k) ^> v(k) still holds; and (iii) relativistic 
scaling region k <C «2, c = v and the Lorentz invariance 
is recovered. 



D. Electromagnetic properties 

Let us discuss the permittivity e(k) and the permeabil- 
ity /j.(k) — 1 + 4ttx (x : magnetic susceptibility). The an- 
alytic solutions obtained from Eq. (fTT]) and /i = l/(ec 2 ) 
are shown in Fig. [3] The momentum scale k can be re- 
garded as the temperature T by T = v(k)k. As noted 
above, the velocity v{k) is the function of the momen- 
tum scale, hence the energy dispersion E{k) = v(k)k is 
a nonlinear function of k. From the definition of k\ and 
K>2, v(ki) = vq and v(K2) — Coo, and the corresponding 
temperatures are estimated as T\ — T{k\) = VqKi and 
T2 = T(n2) = CaoK.2- In Fig. |31 it can be seen that the 
permittivity e(/t) grows logarithmically below T\ while 



4 



the permeability (i(k) decreases below T 2 . The orbital 
magnetic susceptibility x without the electron-electron 
interaction logarithmically diverges as a function of T, 
but in our analysis, the logarithmic divergence is can- 
celed due to the renormalization of v. These contrasting 
behaviors of e and fi facilitate the identification of T\ 
and Ti experimentally. In the zero temperature limit e 
diverges while \i goes to zero, i.e., the perfect diamag- 
netism \ — — l/(47r) is accomplished. 



E. Spectral function 



Im G(to) 




For the physical interpretation of the electron Green's 
function, we should consider the self-energy in Coulomb 
gauge as discussed above. In Coulomb gauged, the elec- 
tron self-energy is 



S(P) = - 
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In principle, ARPES can measure the energy dispersion 
E{k) = v(k — k)k, which shows crossovers at k,\ and k 2 . 
For the spectral function of electrons, the electron field 
renormalization is required, which is given by 



7 2 (w,c, e 2 /e\K) 



1 d\nZ 2 

— K 

2 dn 



e 2 /e 
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From the Callan-Symanzik equation, the electron Green's 
function is modified by the momentum-scale dependent 
functions v(k), a(k), and 72(fc), then we obtain 



G(k,u) 



2 — v 2 {k)k'- 



■ exp 



din 



(18) 

where k! = (uj 1 ', vk') is a four- momentum, A is the energy 
cutoff, and Q is a function determined from a perturba- 
tive renormalization calculation. 

In region (i), 72 = 0, so the Green's function is un- 
changed. In region (ii), k dependence of 72 is rather 
complicated to calculate G(k,u>), so we only consider the 
rclativistic scaling region (iii), where v approaches c and 
the original QED regime is applicable. When we put 
c = v = Coo, 72 (&) is expressed as 



72 (k) 



a{k) 
16tt 2 



Of) 



1 
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on 
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and the perturbative correction for Q is 



g(a(k)) = l + 



a{k) 



FIG. 4. (Color online) Frequency dependence of the density 
of states in region (iii) with Qo = 1. The vertical axis is in 
linear but arbitrary scale. The solid line denotes the result 
with the perturbative correction in Q, while the dashed line 
depicts the result for Q = 1. 



Then, the Green's function becomes 



G(k,u) 



G(a(k)) 



1 



In 



A 2 



3/4- 



12tt 2 

*~ -00 ■- ' -■ 

(21) 

By substituting uj with uj + iO, the imaginary part of the 
Green's function — ImG(fc, w) gives the electron spectral 
function. The electron spectral function has finite value 
for |w| < \k\, while —lmG(k,uj) — outside that region. 
As depicted in Fig. |H the perturbative correction for Q 
gives very small contribution, so we put Q = 1 in the 
following analysis. When the bare coupling constant a 2 , is 
small enough, the spectral function has the approximate 
form 



-ImG(k,u) 



2 c 2 k 2 ) 



327TCoo|fc| \Coolfc 



- LJ / 

(22) 



where the residue a is a constant determined from the 
sum rule. The 5 function peak with finite a means that 
the system remains a Fermi liquid in sharp contrast to the 
(2+l)D case^, while the continuum state for \u\ < c^k 
comes from the interaction as shown in Fig. 2) 



F. Energy gap 

Up to now, we have focused on the critical point 
(to = 0), but the mass to is a relevant parameter. Ex- 
perimentally, the bare mass too can be controlled by the 
concentration x or by pressure P— The RG equation 
for mass to(k) is 



dm(n) 3a(n) 



dn 



8tt 2 



-TO. 



(23) 
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Then, the mass at momentum scale p is 



m(p) = m(A) 



6tt 2 



ln['^ 



9/4 



(24) 



When we neglect the weak singularity with log log mo, 
the solution to Eq. (IM1) is given by m = m [l + 
ln( — )1 9 / 4 , which describes the critical behavior of 

D7T ^ mo /J ' 

the gap as a function of mo oc (a; — x c ) or mo oc (P — P c ) 
with x c (P c ) being the critical concentration (pressure). 



III. DISCUSSIONS AND SUMMARY 

Now we discuss the relevance of the present results 
to the real systems. The velocity vq is estimated at 
vq = 10 6 m/s from the ARPES measurement of the en- 
ergy dispersion^, hence c va cuum/wo — 300. As the di- 
electric constant, we take the typical value £o — 10 2 
for BiSb alloys^. Since c = c va cuum/V£o, c /v Q = 30, 
(co/«o) 2/3 = 10, and a = (c vacuum /w ) x 1/137 = 3 
are obtained. These values give the estimates for K\ = 
10 Ko and Ki being extremely small. Unfortunately, 
it would be difficult to observe the effect of electron- 
electron interaction and the scaling behavior at the ex- 
perimentally accessible temperature in the materials at 
hand. However, there are many candidates for the corre- 
lated topological insulators recently proposed and partly 
synthesized— ~—. The smaller value of vo rapidly (expo- 
nentially) increases K\, which gives the clue to look for 
the appropriate materials to study the scaling behavior 
of the quantum criticality. 

In summary, we have studied the (3+l)D Dirac elec- 
trons coupled to electromagnetic field as the model for 
the quantum critical phenomenon of the transition be- 
tween topological insulator and trivial insulator. The 
RG equations are derived and the two scaling regions are 
identified, (i.e., the nonrclativistic and relativistic scaling 
regions). The Lorentz invariance is recovered in the latter 
case. The physical properties such as the the permittiv- 
ity, the permeability, and the electron spectral function 
have been discussed based on the RG equations. 
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Appendix: Precision of the analytic solutions 

In the appendix, we consider the precision of the ana- 
lytic solutions to the RG equations, especially the valid- 
ity of the relation c 2 v = const. In this section, we define 
x = 1 — v/c, y — c 2 v. The RG equations for x, y and 
e 2 /e are 



dx _ e 2 /e 1 
dn 12-7T 2 y 
dy e 2 /e 



dn 6ir 2 ^ ' 
d(e 2 /e) (e 2 /e) 2 1 



dn 



6ir 2 y ' 
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where 
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From Eqs. (|A~2|) . (jA~3)) . we obtain 

e 2 /e dy 



(A.4) 
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V d(e 2 /e) 
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We assume v < c, i.e. < x < 1, so that < g(x) < 
17- 12 V2 = 0.03. The RHS of Eq. (jA~6|) is small, so the 
relative difference of y is 



(e 2 /e)(K) 



el/eo 



-g( x ) 



(A.7) 



The maximum value g(x) = 0.03 is rarely observed. 
Thus, y — c 2 v can be regarded as a constant to a large 
extent. 
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